Decaying magnetohydrodynamics: effects of initial conditions 
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We study the effects of random homogenous and isotropic initial conditions on decaying Mag- 
netohydrodynamics (MHD). We show that for an initial distribution of velocity and magnetic field 
fluctuations, appropriately defined structure functions decay as power law in time. We also show 
that for a suitable choice of initial cross-correlations between velocity and magnetic fields even order 
structure functions acquire anomalous scaling in time where as scaling exponents of the odd order 
structure functions remain unchanged. We discuss our results in the context of fully developed MHD 
turbulence. 
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I. INTRODUCTION 

In the absence of external forcing the fluid cannot maintain its steady state because of the continous dissipation 
of energy. The decay of kinetic energy in incompressible fluid turbulence has been explored in several studies (see, 
e.g., ||l|-^). These studies, principally numerical, indicate that kinetic energy decays as '^ (i — to)~^ , where to is the 
initial time, i.e. the time when the initial forcing is switched off. The value of t] seems to depend on the initial state. 
Similarly, recent numerical studies of decaying MHD turbulence |^-^ also suggest that the total energy decays as 
power of time. The phenomena of decaying MHD turbulence is observed in many different situations. Star-forming 
clouds are astrophysical examples |g[ . 

A. Decay of MHD turbulence in a shell model 

In Ref. [0 we described a new shell model for MHD turbulence, which has all the conservation laws in the inviscid, 
unforced limit, has the right symmetrires and reduces to the well-known GOY shell model for fluid turbulence in 
absence of any magnetic field. In this paper we use our shell model to study decaying MHD turbulence. We begin 
by solving this set of equations with the external forces and by using the numerical method described there. We 
evolve the shell velocities w„ and magnetic fields 5„ till we achieve the nonequilibrium, statistical steady state whose 
properties are described in Ref. m. Once we attain this state, say at a time t — tss, we use the last set of shell 
velocities and magnetic fields, i.e., Vn{tss) and b„{tss), as the initial condition for our shell-model study of decaying 
turbulence in which we solve our shell-model equations but with all the forcing terms set equal to zero. We then use 
the resulting time series for Vn and bn to determine the total energy E{t) = J^kiWO^^ OP) + (|t*(^' ^)\'^) ^^^ ^^'^ total 
crosshelicity H'^{t) = (|v(k, i).b(k, i)|) where the origin of time is taken to be zero. We find that E{t) ~ H'^{t) ~ t~P, 
with p ~ 1.2; this is illustrated in the log-log plot of Fig.l. As we have noted above, it is relevant for initial conditions 
that have the correct statistical properties of homogenous, isotropic MHD turbulence which are well decribed by our 
shell model. 
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FIG. 1. A log-log plot of decay of total energy and crosshelicity versus time. 



In a recent paper g] the effects of random initial conditions on the decaying solutions of the unforced Navier- 
Stokes equation with hyperviscosity have been discussed: It has been shown that, for an initially Gaussian velocity 
distribution, equal-time two point correlation functions like < u{x,t)u{x,t) > decay as ^ t^^. It has also been 
shown that, for this particular choice of initial conditions, the theory remains finite (in a field-theoretic sense) just by 
renormalisation of the two-point vertex functions. However, it is more relevant to explore the statitical properties of 
decaying fluid turbulence with ordinary viscosity. It is also important to find out how the introduction of magnetic 
fields influences this decay. In this paper we study decaying MHD with ordinary viscosity; we work with both the 
SdMIID equations and a d-dimensional generalisation of the IdMHD model introduced in Ref. M. We show that 
both these models give similar results. In particular, we show that, in both these models, appropriately defined 
structure functions decay with a power law in time; this is modulated by a scaling function whose form depends 
on the initial state. In Section. II, we describe the equations we work with and the intial conditions we use. For 
most of the calculations in this Chapter, we choose initial correlations which are analytic in the low-wavenumber 
limit. We choose ordinary viscosities for both velocity and magnetic fields. We calculate the dynamic exponent z 
and show that it is unaffected by the presence of the nonlinearities, i.e., it remains 2. We next show that if there 
is no initial cross-correlation between velocity and magnetic fields, correlation functions defined at the same spatial 
point decay with power laws in time and, consequently, two point structure functions also decay likewise. Also, 
multipoint correlation functions (and structure functions) exhibit similar behaviour, with simple (and not anomalous) 
scaling in time. In other words, the exponents for high-order structure functions are simply related to those of two- 
point functions. Another interesting feature is that there is no change in the exponents with the introduction of the 
magnetic field (i.e. the exponents are the same for fluid and MHD turbulence); only the amplitudes of the correlation 
or structure functions change. In the Section III we show that in the presence of suitable cross correlations higher 
order correlations (and consequently structure functions) decay more slowly than when there is no cross correlation, 
i.e., they display anomalous scaling. In Section IV we extend and discuss our results in the context of the decay of 
fully developed MHD turbulence. In Section V we conclude and compare results obtained from both SdMHD and our 
simplified model. 

II. DECAYING MHD TURBULENCE: THE CASE OF ZERO INITIAL CROSS-CORRELATION OF THE 

VELOCITY AND MAGNETIC FIELDS 

A. Construction of a simplified model 

The dynamical equations which govern the time evolutions of the velocity and magnetic fields in a magnetised fiuid 
are given by the equations of Magnetohydrodynamics, namely, 

9^ , ^^ Vp (V X b) X b ^2 n /ix 

— + v.Vv = ^ + ^ -^ + i.V2v + fv, 1 

ot p 47r 

dh 

— =Vx(vxb)+AiV^b + fb, (2) 

ot 

where v and b are the velocity and the magnetic fields, i^ and /i arc fiuid and magnetic viscosities, p and p represent 
pressure and density, and fv and fb are external forces. We also impose the incompressibility condition V.v = and 
the divergence-free condition on the magnetic field V.b = 0. In the absence of magnetic fields the MHD equations 
reduce to usual Navier-Stokes equation for fluid turbulence. 

We now briefly present the ideas that go into the construction of the model, which has already been discussed in 
details in ||9[. Since this model is an extension of Burgers equation for fluid dynamics to MHD, it should have the 
following properties: 

• Both the equations should be Galilean invariant. 

• The equations should be invariant under t -^ t; x — > — x, v -^ — v and b ^ b, i.e. the velocity field is odd parity 
and the magnetic field is even parity. 

• The fields are irrotational vector fields in d-dimension ( similar to Burgers velocity field). 

• The equations in the ideal limit i.e. in absence of any viscosity and any external forcing should conserve all 
the scalar conserved quantities of 3dMHD namely the total energy and the cross-helicity. The third conserved 
quantity i.e. the magnetic helicity is zero for irrotational fields . 



The above considerations lead to the foUowing model equations 

dv 1 o 1 T o 

— + -\/v^ + -\/b^^iy^^v + U, (3) 

|^ + V(v.b)=A*V2b + fb, (4) 

whose relation to the 3dMHD equations is the same as the relation of the Burgers equation to the Navier-Stokes 
equation. Here v and b are the velocity and magnetic fields with magnitudes v and b respectively. These fields are 
irrotational, i.e., V x v = and V x b = 0. These equations together conserve the total energy 1/2 /(w^ + 6^), d'^x 
and crosshelicity / v.b d'^x. 

B. Choice of initial conditions 

We have < Ui{k, t — 0) >~< bi{k, t ~ 0) >— by homogeneity and isotropy. We choose the initial correlations to 
have the simple analytic forms 

<w,(fc,t = 0)wj(fc',t = 0) >^Duuhkjd{k + k'), (5) 

< b^{k,t = 0)bj{k\t = 0) >= Dbbhkj6{k + k'), (6) 

and 

< v^ik, t = 0)^bJ{k\t = 0) >= D^j{k)6{k + k') (7) 

In fc-space the SdMHD equations are 

-i^p\{p) + Mp) (8) 

dMp) = -{ip X Y^ v{q) X 6(p - q))^ - ^ip\{p) + ^^(p) (9) 

The quantities My;(k) — Pijki + Pukj appear to include the incompressibility, Pij{k) = (Sij p^) is the transverse 

projection operator and wavevector arguments p, q, k indicate spatial Fourier transforms. The initial conditions we 
use for SdMHD are 

<Ui{k,t^O)uj{k',t^O)> ^ Duuk^Pzj{k)S{k + k'), (10) 

<bi{k,t = 0)bj{k',t^O)> ^ Dbbk^P,j{k)6{k + k'), (11) 

<Uiik,t = 0)bj{k',t=^0)> = D,j{k)S{k + k') (12) 

which ensure the divergence- free conditons, namely, 

h < u,{k)uj{-k) > =0 (13) 

fc, < 6,(k)6,(-k) > ==0. (14) 

We discuss the structure of Dij later. For the time being, we choose Dij ~ 0. These simple initial conditions, which 
are analytic in the fc ^ limit, arc enough to elcidate the issues raised in this Chapter. 

C. Calculation of effective correlation and structure functions within a one-loop perturbation theory: 

temporal scaling 

The formal solutions of the Eqs. (3) and (4) are given by 



ui{k,t) = e-''''^'ui{k,0) - ^ki I dt' f d'^pexp[-i^{t - t')k'^]exp{~i^t'[p'^ + {k ~ p)^]} 



'0 

„2 I /u „\2l 



[um{p, 0)u„(k - p, 0) + cxp{-fit'[p^ + {k- pY]} X 

6m(p,0)6™(k-p,0)] + ... (15) 

6/(k,i) = e"'''^'*6,(k,0)-ifc, / dt' f d'^pexp[^fi{t - t')k^]cxp{-t'[i^p^ + n{k - pf]} 

K(p,0)6,„(k-p,0)] + ... (16) 

where ellipsis refer to higher order terms Similarly the formal solutions of the 3dMHD equations (j^ and (|9|) are given 
by: 

ui{k,t) = e-^'='*uKk,0) - ^Mij, f dt' j d'^pexj>[~v{t - t')k^]exp{-vt'[p^ + {k ~ pf]} 

[uj(p,0)u5(k-p,0) 

-exp{-A*t'[p2 + (A;-p)2]}6^.(p,0)6,(k-p,0)] + ..., (17) 

6,(k,t) = e~'^'='*5/(k,0)-iA:,„ j dt' f d'^pexpl-fiit - t')k^]exp{-t'[:yp^ + fi{k ~ pf]} 

[u„(p, 0)6,(k - p, 0) - ui{p, 0)b,n{k - p, 0)] + ..., (18) 

where the ellipsis in refer to higher order terms in the series. We use bare perturbation theory, i.e., after truncating 
the series at a particular order, all v and b are to be replaced by their bare values, i.e., by G"u(i = 0) and G^J){t = 0), 
respectively, where G'"(k, i) = e^"^ * and G\{k,t) — e~^^ *. This is similar to Ref. |lO). In the steady state of fluid 
or MHD turbulence, the dynamic exponent z is determined by the form of the correlations of the driving noise and 
dimensionality of the space. For example, when driven by noise correlations of the form k^^ (in 3dNS), one finds z 
(within a 1-loop dymanic renormalisation group calculation) to be 2/3 [O. What is z when we study decaying MHD 
turbulence? We calculate this within a 1-loop diagramatic perturbation theory: In the bare theory (i.e., when the 
nonlinear terms are absent) z = 2. If the one-loop corrections to the response functions (i.e., to the viscosities) have 
infrared divergence at zero external frequency, then z < 2; otherwise z = 2. The diagrams shown in Fig. 2 contribute 
at the one-loop level in renormalisation of the fluid viscosity v, 

D„„ 





FIG. 2. 1-loop diagrams contributing to the renormalisation of v. 
Similar diagrams renormalise /i. The effective 1-loop fluid viscosity is given by 

1 1 , Dun /■ ,3 PiPmk.(k-p) 1 fc2 



vk'^ vk'^ vk"^ J 2i^k.p 2vp'^ — 2i'k.p i>{p^ + (k — p)^ ' 

1 .-, k"^ M , Dbb /■ 3 pip.mk.{k-p) 1 vk"^ 

Kl - ^^71:2 )] + -T2 / rf P TTi— [^7-2(1 



ivp"^ 2vp'^ vk"^ J 2/L/,k.p 2/ip-^ 2/xp^ 

1 , l^k^ M . N 

r(l-^^^)]- (19) 



2^p2 ^ 2^p2 

We see that, in the long- wavelength limit, 1-loop corrections have no infrared divergence at d = 3; in other words z 
remains 2. In the case of the 3dMHD equations, the forms of the integrals remain unchanged, only changes appear in 
the coefficients of the integrals where, in place oipiPm, appropriate projection operators are present. Thus again there 



is no diagramatic correction to viscosities in the long-wavelenght limit in 3dMHD. Hence the exponent z is unaffected 
by the nonhnear terms. Therefore, without any loss of generality, we can put v = \x = 1. 

We now calculate the effective D^u and D\,\y for d = 3 at the one-loop level. The diagrams shown in Fig. 3 contribute 
at the 1-loop level to -Duu, These contributions are proportional to 
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FIG. 3. 1-loop diagrams contributing to the renormalisation of D^u and D^f, respectively. 

hk, {Dl^ + Dl) J d3pPW_S ^ k,k, [Dl^ + Dl) f £p^ (20) 

which does not have any infrared divergence in any dimension. The one-loop diagram contributing to Dbb has the 
same form. For 3dMSD the forms of the one-loop integrals do not change; only the coefficients have appropriate 
factors of projection operators instead of kikj. On purely dimensional grounds for d = 3, 

(w,(x,t)wj(x,t)) - f d^k{u^{k,t)uj{-k,t)) - Dun f d^ke-"''' *€-"''' *hkj ~ Duut^S^j (21) 

and 

(Wx,i)Mx,«)) ~ /.»*(»*, <)M-k..»a. -^ /.=*»-%-■=■«, ^ a.,^.«, ,22, 

Similar quantities obtained from SdMHD also show the same exponents but with different amplitudes f|. Note that 
the temporal dependences can be scaled out of the momentum integrals easily as the lower and upper limits of the 
momentum integrals can be extended to and oo, respectively. 

It is easy to see that all multi-point correlation functions also do not exhibit any infrared divergence in three 
dimensions: Consider, e.g., 3-point functions. Even though (u(ki, t)Mj(k2, t)ws(k3, t)) = at t = 0, it becomes non- 
zero at i > as can be easily seen at the one-loop level. (u(ki, t)wj(k2, t)Ms(k3, t)) has a tree-level contribution and 
next a one-loop contribution. The tree-level part being proportional to (ui(ki,0)wj(k2, 0)us(k3, 0)) is zero because 
we have chosen a Gaussian distribution for the initial conditions. The one-loop diagrams are shown in Fig.(H). 



^This is easy to understand: Calculation of the exponents through the one-loop method depends upon the power-counting 
of the one-loop integrals which are same for both the 3dMHD equations and our simplified model as the projection operator 
Pij that appears in the 3dMHD equations is dimensionless and the nonlinearities in both the equations do not renormalise due 
to the Galilean invariance of our simplified model and 3dMHD equations (which in turn keeps the power-counting structure 
same). Amplitudes of course depend upon the specific models. Similarly, energy spectra obtained in the steady-state from 
diff'erent models which give same power-counting in the one-loop integrals, when driven by stochastic forces with correlations 
having same scaling, have same scaling propertes, see, for e.g. Refs. [p|,[Ll[. 
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FIG. 4. One-loop diagrams contributing to {u(ki, t)uj(k2, t)us(k3, t)). 
The 1-loop contribution is however nonzero. It has the form 

{u(ki,t)uj{'k2,t)us(k3,t)) r^ exp[iy{-kl - k^ - (ki + 'k2)^)t]kiik2ji-kis - fc2s)Am«, 
where Duuu is the value of the one-loop integral in the limit /ci, fc2 — *■ 0, i.e., 



D„ 



d'^q{D: 



uu ' " ^bb) 



(23) 



(24) 



which is infrared (IR) finite. We can interpret this one- loop result by saying that the nonlinearities in Eqs.(pf) and (H) 
are absent but there are nonzero initial three-point correlations at t = 0, i.e., the effective linear system behaves as if 



Hence 



< Ui(ki, 0)uj(k2, 0)us(k3, 0) > = kiik2jk3sDuuu5{'ki + k2 + ka), 
< 6j(ki, 0)6j(k2, 0)6^(k3, 0) > = kuk2jk3sDbbbS{ki + ka + kg). 



< Ui{-Xi,t)Uj{-X2,t)Us{-X3,t) > 

= / d^fcid^fc2exp[iki.(xi - X2) +k2.(xi - X3)](ui(ki, t)Mj(k2, t)ws(-ki - k2,i) > 
= / d^kid^k2exp{-kl - kj - (ki + k2)^)iexp[iki.(xi - X2) + k2.(xi - X3)] 

,Xi-X2 X1-X3, 



(25) 
(26) 



kiik2jk3sDu 



^ -^ijsx' 



Vt ' Vt 



-)• 



(27) 



The last step can be easily obtained by using the dimensionless variables j/i = kft, 1/2 — A:|i. 

Next we consider the four-point functions: We concentrate on (^^(xi, i)uj(x2, t)us{:ii.3, f)ii;(x4, t)); the diagramatic 
representations up to the one- loop contributions are given in Fig.(p|). 
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FIG. 5. Diagramatic contributions to (ui(ki,i)itj(k2, i)us(k3,i)ui(k4, t)) upto one-loop level. Diagram (a) is the bare 
four-point function due to the initial two-point correlations, diagrams (b) and (c) are one-loop contributions. 



The one-loop integrals in the long- wavelength limit, are 

J Jy fJ, 



(28) 



These one-loop contributions can be interpreted as in the previous case of three-point functions i.e. as if there were 
initial four-point correlations in absence of the non-linearities: 

(Ui(ki,0)wj(k2,0)u5(k3,0)u;(k4,0)) ee DunuukiihihshiSCki + k2 + ka +k4). (29) 

Similar results can be easily obtained for other four-point functions involving velocity and magnetic fields. Hence, 

< Ui(xi,t)Mj(x2,t)-Us(x3,i)u;(x4,t) > 

~ Dl^ / (fki(fik2kiikijk2sk2i exp[-2z^(fcj -|- fc|)t] exp[iki.(xi - X2) + k2.(x3 - X4)] 
+ Duuuu / (fki(fk2d^k3kuk2jk34~kii ~ k2i - kai) exp[-i^{fcj -f fc^ + fc| + {kf + kl + kl)}t] 



(30) 



exp[iki.(xi - X2) + k2(x2 - X3) + k3(x3 - X4) 



n2 t-^ +n t-13/2 

uu "^ ^-^uuuu^ 



In the infinite-Reynolds-number limit, i.e., in the limit i^, /i — > 0, D^uuu -^ 00 and hence the second term in the last 
line of Eq.(pG) dominates. As in the case for the three-point functions four-point functions, are of the form 



< 



Ui{xi,t)Uj{x2,t)Us{x3,t)ui{x4,t) >r^t ^^^'^A^jsl{ 



Xi - X2 X2 - X3 X3 - Xi 



Vt ' Vt 



Vt 



(31) 



Here Aijg and Aijsi are third and fourth rank tensors respectively, and < Ui(x, t)** >~ t i3/2_ rp^ study decaying 
MHD, we extend the definiton of equal-time structure functions as follows: 



and 



5;'(r,i)=<|[u,(xi,t)-ii,(x2,t)]r> 
5,"(r,t)=<|[&,(xi,i)-6,(x2,i)]|">, 



(32) 
(33) 



7 



where r = |xi — X2I and t denotes the time that has elapsed after the forcing terms are switched off. We now look at 
the general forms of these structure functions. Consider the three-point function first: 

^3(r,t) = (|[u,(xi,t)-«,(x2,t)]H. (34) 

On expanding the expressions for, S!^{r,t) we get the following kind of terms: (i)('Ui(xi, t)'^) and (ii)(Mi(xi, i)^'Ui(x2, i)). 
The first one is given by 

< u,(xi,t)3 > = / (fki(fk2kuk2i{-ku - fc2j)exp[-{fc2 - fc^ - (ki +\^2f}t]Duuu 

^ t-9/^ (35) 

whereas the second term has the form 

< 'Uj(xi, t)2u,(x2, t) > = / (fkid^k2 kuk2t{-ku - k2^) cxp[-{fc? -kj- (ki + ka)^}^] 

exp[iki.(xi - ^2)]Duuu - t-'^'^4){r/Vt), (36) 

with r = |xi — X2I. Thus, it is easy to see that S'^{r,t) has the structure 

SUr,t)^t-'/'F^{r/Vt), (37) 

with F^{r/^/t) is a function with dimcnsionless argument r/^/t. We can easily extend our calculation above to see 
that 

^:(r,t)^t-^"^«(^), (38) 

5,"(r,i)^t-^"^^(^), (39) 

where F^ and Fj^ are functions with dimcnsionless argument r j \/t. It is obvious that all these higher-order corre- 
lation functions do not exhibit any anomalous scaling in time; their temporal behaviours can be obtained by simple 
dimensional arguments from the behaviours of the two-point functions. It is easy to see that 

gn ^ i-(n-l)3/2-„/2^.(_l_)^ (40) 

^„^^-(„-l)3/2-„/2^.(_^)^ (41) 

i.e., z^ ~ z^-^ = (n — 1)3/2 + n/2, at the level of our one- loop approximation. 

III. NON-ZERO INITIAL CROSS-CORRELATIONS 

A. Form of the initial cross correlations 

What happens if we make Dij non-zero? Let us first establish the form of Dij-. Recall that < Ui(xi, i)6j(x2, i) > 
has odd parity (since u is a polar vector whereas b is an axial vector). Hence it is purely imaginary and is an odd 
function of k. (From now on we will work only in three dimension.) Thus we choose 

where e^p is the totally antisymmetric tensor. This is the simplest analytic form with the desired structure. We 
consider the cases a — and 2. 



B. The temporal decay of correlation and structure functions 



It is easy to see that the value of z docs not change if D^h 7^ as no new diagram appears. However, an additional 
diagram (Fig. 4. 6) contributes to the effective D^u and Dbb if D„(, 7^ 0. 




FIG. 6. 1-loop diagram contributing to the renormalisation of Duu and Dhh, which arises when D^h 7^ (see text). 
This diagram yields the following contribution 

Sq 26snq_s<ln<f 



kiik 



2i 



16fc4 



D^ 



(27r)3 



where we have used the identity 



^Ims^lmn \^ ^)^s 



2S, 



(43) 



(44) 



Naive power counting shows that this contribution is finite for d = 3, hence Duu and Di,h are not renormalised even 
if Dub 7^ 0, for both a = and 2. The same conclusion holds for SdMHD as the one-loop integral remains finite; 
only the coefficients change because of factors of Pij{k). It is quite interesting to note that for, Dub 7^ 0, there is a 
diagramatic correction to Dub', 

D„„ 





the form of the correction is 



FIG. 7. 1-Ioop diagramatic correction to Dub 



kiik2j j^ j^ f d^q <f>^hnsqsqiqn 



IQk^ 



(27r)a 



(45) 



however this correction vanishes for both a — Q and 2 since the integrand is odd in q. Hence there is no renormalisation 
of Dub and two-point structure functions Eqs.(4.22) and (4.23) scale as t"^/'^ just as they do when Dub = 0. 



C. Higher-order correlation and structure functions: anomalous temporal scaling 

We now consider the ri-point correlation functions with n > 2. None of the three-point functions have any infrared 
divergence. Hence they exhibit simple scaling, i.e., < Ui{x,t)^ >~ i^^/^. Let us consider the four-point functions 
namely, 

< Mi(xi,i)Uj(x2,i)'Us(x3,t)Mi(x4,i) >, 

< bi{xi,t)bj{x2,t)bs{-K3,t)bl{x4,t) >, 

< Ui{xi,t)Uj{x2,t)bs{x3,t)bl{x4,t) >, 

< Ui(xi,i)Uj(x2,i)Us(x3,t)6i(x4,t) > . 

(46) 

At the 1-loop level, all of them they will have finite parts (coming from initial auto correlations of velocity and 
magnetic fields). However, the first three of them will have new, infrared-divergent, one- loop contributions if D^t 7^ 0. 




oD 



ub 



oD 




ub 




Us Uj bs ^\ ^s ^i 

FIG. 8. Divergent one-loop contributions to (Mi(xi, t)uj(x2, t)us(x3, t)u;(x4,i)), (ui(xi, t)uj(x2,i)6s(x3, t)bi(x4, t)), 
{6i(xi,i)&j(x2, t)fes(x3,i)6i(x4,i)) respectively (see text). 

The analytic structure of the divergent 1-loop integral is illustrated for one of these correlation functions in Fig. (0) 
whose contribution is 



(ui(ki)uj(k2)us(k3)u/(k4)) 



/c8 



d^P PmPnPfPtP° 

/rt \o^acm^cdn^def^eat o 

(27r)3 p8 



Isotropy yields 

PmPnPfPt = -^P'^[SmnSft + SmfSnt + SmfSfn] 

SO summing over (m, n) and (/, t) separately we obtain 

1 



A 



1/15 



d2+2d 

for d = 3. Hence naive power counting for the integral in Eq.(|47|) yields, for a — Q, 

kuk2jks,skii ^4 [ d^pp'^ 



15/fc8 '^^b 



p° 



kuk2jk3sk4iA 



(47) 



(48) 



(49) 



(50) 



where A is a momentum scale arising from the lower limit of the integral. Hence it is linearly infrared (IR) divergent 
for a = 0. For a = 2 there is no IR divergence as is the case with D^b = 0. We analyse the case a = in detail: 
There is a disconnected diagram also that contributes to {ui{xi,t)uj{x2,t)us{x3,t)ui{x4,t)), which yields 
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A, 



Dt 



(51) 



3 ' "'4J1 



(52) 



This is finite in the fc ^ hmit, so the contribution from the connected part will dominate in the /c — > Hmit. If 
we define the effective four-point correlation to be ex Duuuu ~ Duuuu^~^ (with the choice k = yHt 
where -D^uum i^ J^^^ a constant) then 

(ui(ki,0)Wj(k2,0)Ms(k3,0)w;(k4)) ^ fciifc2ifc3sfc4ifc~^£)°„„^5(ki +k2 +k3 +k4). 

Hence one obtains 

< -u,(x, tf > r^ j cfqid^q2(fq3exp{-kl - k^ - k^ - (fci + fca + fca)^))^ 
xkiik2ik3i{kii + fc2, + hi)D'^uuuk~^ 

.-3d/2-2.-l/2 +-1/2 ^ „ (^ ,\4 ^ l 
~t t ' ^t ' <Ui[X,t) > \D^t=0- 

Similarly (6.(x,i)4) ^ t-i/2(„,^(x, t)4)|^^^^Q^ < u,(x, t)26,(x, f)^ >^ i-1/2 < m,(x, t)26,(x, t)^ > | 



(53) 

_D„[,=o- The other 

four-point functions are (dimensionally) linearly divergent at the 1-loop level; however, the 1-loop integral vanishes 
becauses the integrand is odd in p. Hence they will scale as t"^^^"^. It is now obvious that 



< K(xi,t)-^i,(x2,i)]'> ^ i-l/V13/4A„( 



Xl - X2 - 

tl/2 ' 



- r 1/2 < [u,{xi,t) - M,(X2, t)]^ > b„ 

< [6.(xi, t) - &.(X2, i)]' > - t-^/'t-^'^^A,{^^ ""^ ^ 



tl/2 



•-) 



-1/2 



< [ui(xi, t) - Ui{:x.2,t)]^ > |d„,=o- 



(54) 



(55) 



Hence we see that these structure functions exhibit anomalous dynamical scaling if Dub 7^ 0. Similarly, it follows that 
the even-order exponents are 



L = ^2„ = ;^(2n-1) + : 



'2n ' 



n> 2, 



where zf"° = (2n — 3)/2 is the anomalous part of the exponent Z2n- By contrast the odd order exponents 

z^ = zf^ = (n - 1)3/2 + n/2, for any odd n, 
do not have an anomalous part. In Fig.(ph we present plots of z„ versus n for both even and odd n. 



(56) 
(57) 



N 




n 

FIG. 9. A plot showing the dependence of Zn over n; the line gives the exponents for the odd order structure functions for 
all values of D^^t and for the even order structure functions when D^t = 0; we also show 24 and zq in this plot (circles) when 
£>„(,/ (for a = 0). 
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IV. EFFECTS OF INITIAL CONDITIONS WITH SINGULAR AMPLITUDES FOR THE VARIANCE 

In this Section we choose initial conditions of the following form: 

<u,{-k,0)u,{k',0)> =Duu^5{k + k'), (58) 

<fe,(k,0)6j(k',0)> ^Dbb^S{k + k'), (59) 

and 
<u,(k,0)6,(k',0)> =iD„be,,pA_j(k + k'). (60) 

We chooseQs > y. The effective fluid viscosity i), and two-point correlations Duu and D^u are given by 



i>k^ vk"^ vk? J v[p'^ + fc2 — (fc — p)2] vIp^ + (fc — p)^ — k'^] p^ + (/s — p)'^ ' 

1/1 ^%, n ^ , 1 /■ d^PP^'PtPjDbb . 1 

r(l ^)k.(k-p) ' ' ^ r 



2t'p^ p'^ vk"^ J /i[p^ + fc^ — (fc — p)'^] pl[p^ + (fc — p)^ — k"^] 

p 1 fc2 
(1 9 n ^)~^^^i^ ^)k.(k-p), (61) 

^™ = ^n« + A^^, (62) 

Dbb = Dbb + A-^, (63) 

Dub = Dub, (64) 

where A is just a numerical constant. The self-consistent solution is given by :^, /i ~ k~^ , with the long-wavelength 
behaviours of Z)„„, Dt,f, and D^b unchanged. Now we can easily calculate the decay of the total velocity and magnetic 
energy, and cross helicity; these are: 

6-23 



E''{t)r^t-^^, (65) 

E^{t)^r^^, (66) 

and 

H%t)^t-^^. (67) 

Here E" = Jd'^kl/2 < \v{k)^ >, E'' = Jd'^kl/2 < \b{ky >, H^ = / d'^fc |v(k).b(k)| are total velocity-field and 
magnetic-field energies, and total crosshelicity respectively. While calculating the temporal decay of these quantities 
it has been implicitly assumed that the lower and the upper limits of the momentum integrals can be extended upto 
and cxD, respectively; if this not possible temporal behaviour cannot be obtained by scaling (see next Section). 

A. Decay of fully developed MHD turbulence 

In fully developed SdMHD turbulence one finds, in the inertial range, 

{v.,{k,t)vj{k',t)) ^ P,,(k)fc-"/35(k + k'), (68) 

{b,{k,t)bj{k',t)) - P,j{k)k-^^/^S{k + k'). (69) 

Of course, a study of higher-order correlation functions reveals that the velocity and magnetic- field distributions are 
not Gaussian ||7|. If, however, for the purpose of analytical tractability, we assume that the initial distributions are 
Gaussian with the variances as (for our Burges-type model) 



li y > s then one- loop effective D^u and Dtb will scale as A; ^. 
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(w,(k,t = 0)uj{k',t^ 0)) = Amfc^fc,|A:|-"/'^^ (70) 

(&,(k,t = 0)6,(k',t = 0)) = A6fc,fc,|fcr"/3-^ (71) 

such that I < u{k,t = 0)u{-k,t = 0) > | -- | < b{k,t = 0)b{~k,t = 0) > | -- fc~"/^ as observed in fully developed 
MHD turbulence, we can calculate z and effective two-point correlations as we did above in Sec (4.2.3). For simplicity 
we give details only for our simple Burgers- type model [Eqs.(^||)]. The effective v, D„„ and Dtb are given by 

1 _ 1 , Duu f,d PlPmp-^^/'-^M^-p) , 1 ,, k^ 



vk"^ vk'^ i/k^ J 2i^k.p 2i/p'^ — 2t/k.p iy{p^ -|- (k — p) 



2iyp'^ 2vp'^ uk^ J 2/xk.p 2^p'^ 2/ip2 

1 , vk"^ ,, , , 

-•(^-^TTZ^)] (72) 



2/xp2 2/ip2 

^«" = A» + ^^+^;^ (73) 

Here (5 is a numerical constant. The self-consistent solutions of these equations are given by i^, /i '-^ k~^/^ and 
Dun = -Dftf,, implying that z = 2/3; this is the same as in the steady state. With this self-consistent solution, 
we see that the 1-loop corrections to Z)„„ and D^fj are as singular as the bare quantities. Hence we ignore their 
renormalisations. When I?„b is non-zero, it will contribute to the 1-loop corrections of !)„„ and D^b with the same 
low-wavcnumbcr singularity. Dub itself does not have any diagramatic correction, as in the case discussed in Section 
III. 

The decay of two-point correlation functions will be of the form 

(^.(x,^)^)- /rf^fc ea;p(-2i/fc2/3i)fc-22/3 (75) 



It is easy to see that, because of the infrared divergence of the integral, the lower limit cannot be extended to 0, 
so it will depend upon a lower cut-off (the inverse of which is the integral scale). Consequently the i-dependence 
cannot be scaled out as in the previous case. This problem is reminiscent of the sweeping divergence that appears in 
a perturbative DRG calculation of fully developed turbulence or MHD turbulence |^,|2| . In fact this problem begins 
when < Mi(k, 0)uj(k',0) >~ fc"*""""^ and y = 3; then the above integral becomes log-divergent. This type of initial 
state can be prepared by randomly stirring the system with a stochastic force with a variance k'^^y^'^ with y = 3. 
In a DRG analysis, sweeping divergences appear when y > 3 (but K41 scaling obtains for y = 4) |l2]. Higher-order 
correlation functions have more severe divergences, hence their decay also depend strongly upon the integral scale. 

V. CONCLUDING REMARKS 

We have examined the decay of the total kinetic and magnetic energies, and the crosshelicity in our shell model. 
These quantities decay with power laws in time. In order to understand this behaviour analytically, we have employed 
a one-loop perturbation theory. In the first part of our calculations, we have chosen random initial conditions with 
the statistical properties given in Eqs.(5),(6) and (10), (11) (with no crosscorrelations). We have shown that such 
initial conditions lead to power-law decay of temporal correlations and appropriately defined structure functions. For 
initial correlations analytic in the long-wavelength limit no higher order structure functions show any anomalous 
scaling. However, on introducing crosscorrelations < Ui(k, 0)&j(k', 0) > (Eq.4.43) higher-order structure functions 
exhibit anomalous scaling. This highlights the importance of crosscorrreations in the statistical properties of MHD 
turbulence. We have also extended our calculations to the case of initial correlations with singular variances and 
discussed the difficulties in the decay of fully developed MHD turbulence. We have argued, that these difficulties are 
related to the sweeping effect which arises in the study of the stochsatically driven MHD turbulence (in the statistically 
steady state). It is also worth noting that we have obtained similar temporal behaviour from both SdMHD and our 
simplified Burgers-type model. It would be interesting to see if our results can be checked experimentally and 
numerically. 
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